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Abstract. We give the basic definition of algebraic entropy for lattice equations. 
The entropy is a canonical measure of the complexity of the dynamics they define. 
Its vanishing is a signal of integrability, and can be used as a powerful integrability 
detector. It is also conjectured to take remarkable values (algebraic integers). 

The analysis of discrete dynamical systems, in particular the measure of their 
complexity, and possibly the detection of their integrability is a huge subjectf, 
originating in the work of Poincare. It contains the study of the dynamics of rational 
maps, already a vast topic of research. It also contains the study of lattice equations, 
which are to maps what partial differential equations are to ordinary differential 
equations. Our purpose here is to extend the notion of algebraic entropy, already 
widely used for maps |2] [HI Qj, and recognized as an unmatched integrability 
detector j^, to lattice equations as in [S], thus introducing a measure of complexity 
for higher dimensional discrete dynamics. 

We briefly describe the setting, the space of initial data, the evolutions, and 
deflne the related entropies. We give examples, some integrable, some not integrable, 
showing how to extract information about the global (and asymptotic) behaviour of 
the system from a few iterates. We formalize and confirm the results of (HI- We also 
present some conjectures and perspectives. 

1. The setting 

Consider a cubic lattice of dimension D. The vertices of the lattice are labeled by D 
relative integers [ni, 712, . . . , n^i]. To each vertex is associated the variable y[ni,n2,...,nD]- 
We are given a defining relation, which links the values of y at each corner of the 
elementary cells (square for D = 2, cube for D = 3 and so on). 

We will suppose that the defining relations allow to calculate any corner values 
on a cell from the 2^ — 1 remaining ones, and that the value is given by a rational 
expression. This implies that our defining relations are multilinear, which covers a 
large number of interesting cases (see for example [3IH] and more in the last sections). 
This restriction may be partially lifted. 

As an illustration, consider the two-dimensional case of the square plane lattice. 
The elementary cell is a plaquette shown in Figure (1). 



t We will not dwell here upon the definition of integrability 
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Figure 1: elementary cell in two dimensions 
The defining relation is a constraint of tlie form 

/(y[rii,n2]; y[iii + l,n2]i y[rii,n2 + l]' y[ni + l,n2 + l]) = 0' 

We will use some specific examples later. 



(1) 



2. Initial conditions 

For the sake of clarity, we will concentrate on the D = 2 case, but everything we will 
say generalizes straightforwardly to higher dimensions. 

In order to define an evolution we have to specify initial conditions. From the 
form of the defining relation, it appears that the values of y have to be given on some 
"diagonal" of the lattice. 

The space of initial conditions is infinite dimensional. 

When D — 2, the diagonals need to go from [ni = — cx3,n2 = — oo] to 
[tli = oo,7i2 = oo], or from [ni = —00,712 = +00] to [771 = 00,772 = —00]. We 
will restrict ourselves to regular diagonals which are staircases with steps of constant 
horizontal length, and constant height. Figure (2) shows four diagonals. The ones 
labeled (1) and (2) are regular. The one labeled (3) would be acceptable, but we will 
not consider such diagonals. Line (4) is excluded since it may lead to incompatibilities. 
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Figure 2: diagonals 
3. Restricted initial conditions 



Here again we use the D = 2 example, to make things simple. Given a line of initial 
conditions, it is possible to calculate the values y all over the -D-dimensional space. 
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We have a well defined evolution, since we restrict ourselves to regular diagonals. 
Moreover, and this is a crucial point, if we want to evaluate the transformation formula 
for a finite number of iterations, we only need a diagonal of initial conditions with finite 
extent. 

For any positive integer N, and each pair of relative integers [Ai, A2], we denote by 
^[Ai''a2] ' ^ regular diagonal consisting of N steps, each having horizontal size h = | Ai | , 
height h = IA2I, and going in the direction of positive (resp. negative) nfe, if Afc > 
(resp. Xk < 0), for k = I..D = 2. See Figure (3). 




Figure 3: various choices of restricted initial conditions 

Suppose we fix the initial conditions on aI^^^'^^^^j. We may calculate y over a 
rectangle of size {Nli + 1) x + !)• The diagonal cuts the rectangle in two halves. 
One of them uses all initial values, and we will calculate the evolution only on that 
part. See Figure(4). 



(3) 

Figure 4: range of initial conditions 2 



4. Fundamental entropies 



We are now in position to calculate "iterates" of the evolution. Choose some restricted 
diagonal A|_^''^^j. The total number of initial points is g = N{li + I2) + 1. 
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For such restricted initial data, the natural space where the evolution acts is the 
projective space Pq of dimension q. We may calculate the iterates and fill Figure (4), 
considering the q initial values as inhomogeneous coordinates of Pq . 

Evaluating the degrees of the successive iterates, we will produce double sequences 
of degrees. 

The simplest possible choice is to apply this construction to the restricted 
diagonals Aj^j' which we will denote A^'', A^|^\ . . . , and call them fundamental 
diagonals (the upper index (N) is omitted for infinite lines). 

The pattern of degrees is then of the form 
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(2) 



To each choice of indices [±1,±1] we associate a sequence of degrees d±j.. 
Definition: The fundamental entropies of the lattice equation are given by 

e±± = lim — loq{d'+l). (3) 

n^oo n 

Claim: These entropies always exist [21, because of the subadditivity property 
of the logarithm of the degree of composed maps. 

The fundamental entropies correspond to initial data given on diagonals with 
slope +1 or —1, and evolutions towards the four corners of the lattice, as shown in 
Figure (5). 
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Figure 5: fundamental evolutions on a square lattice 

These four entropies do not have to be identical (see section 
When the entropy vanishes, the growth of the degree is polynomial, and the degree 
of that polynomial is a secondary characterization of the complexity. 
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5. Subsidiary entropies 

We may also define entropies for the other regular diagonals. They correspond to 
initial data given on a line with a slope different from ±1. They are useful in view of 
the various finite dimensional reductions presented for example in P]. 



As an example, the pattern of degrees for A 

1 



(3) 



-1,21 



looks like 



d(3)[l] = d(6)[2] 

(i(5)[2] 
d(3)[2] 

(i(2)[2] 

d(i)[2] 
1 



(4) 



The sequences of degrees we will retain are the border sequences {d^") 
p = 1,2, seen on the edges of the domain. There are as many such border sequences 
as there are dimensions in the lattice (here D=2). The index in bracket refers to 
the direction of the edge considered. This leads to subsidiary entropies e^Xi Ad][H) 
u = l...D: 



1 



e[Ai,....AB]H = lim - log(drr] Ax,iM) 



(5) 



6. Explicit calculation 

A full calculation of iterates is usually beyond reach. We can however get explicit 
sequences of degrees by considering the images of a generic projective line in Pq, as 
was introduced in Q], making a link with the geometrical picture of 0: 

Suppose we start from a restricted diagonal A[f^\ , . It contains q — N{li+l2) + 1 



vertices Vi, 
form: 



[Ai,A2 

Vq. For each of these q vertices, we assign to y an initial value of the 



oik 



1 



(6) 



where uq, /Sq and ak, Pk, {k = l-.f?) are arbitrary constants, and x is some unknown. 
We then calculate the values of y at the vertices which are within the range of Aj-^-'^_^j . 
These values are rational fractions of x, whose numerator and denominator are of the 
same degree, and that is the degree we are looking for. 

The next step is then to evaluate the growth of degrees. One very fruitful method 
is to introduce the generating function of the sequence of degrees 



oo 



k=0 



(7) 



and try to fit it with a rational fraction. 

The remarkable fact is that it again works surprisingly well, as it did for maps, 
although we know that it may not always be the case ^OI- This means that we can 
often extract the asymptotic behaviour measured by ^ and © just by looking at a 
finite part of the sequence of degrees. 
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The existence of a rational generating function with integer coefficients for the 
sequence of degrees implies that it verifies a finite recurrence relation. For maps this 
can sometimes be proved through a singularity analysis. A similar singularity analysis 
should be done here, but it is beyond the scope of this letter. 

We have explored a number of examples. We found examples with non zero 
entropy, and examples with vanishing entropy either with linear growth either with 
quadratic growth, up to now. 



7. Example 1: the deformed cross-ratio relation 

Take as a defining relation: 

„ _ {y[ni,n2] - 0'y[ni + l,n2]) {y[ni,n2 + l] - by[ni + l,n2 + l]) _n I o\ 

Jdcr - 7 TT -: T - S - U (8j 

l,y[ni,n2] ~ Cy[ni,n2 + 1]) l,2/[ni + l,n2] ~ " 2/[ni + l,n2 + l] j 

This relation is based on a deformed version of the cross ratio of the four corner values. 
It is known to define an integrable lattice equation for a = b = c = d ^l] . 

At generic values of the parameters we get the following explicit values for c?^^ 
(the defining relation being very symmetric, the four sequences are the identical): 



= {1, 2, 4, 9, 21, 50, 120, 289, . . .}. (9) 
This sequence is fitted by the generating function 

The entropy is the logarithm of the inverse of the modulus of the smallest pole of ^(.s) 



e 



dcr 
±± 



log(l + V2). (11) 

We have calculated a number of subsidiary entropies for various values of [Ai,A2]. 
They all give sequences which can be fitted with rational generating functions. The 
entropies depend on the "slope" <t = A2/A1: large a give larger e[Ai,A2][l] ^-^d smaller 
e[Ai,A2][2], and conversely. As an example, we see that e[p,i][l] is the inverse of the 
logarithm of the smallest modulus of the roots of 1 — s — — s^^^ . There is an 
interesting interplay between the various fundamental and subsidiary entropies. 

For the known integrable case (parameters a = b = c = d) ^I], we find the same 
sequence for all four A±±: 

{4"±} = {1> 2, 4, 7, 11, 16, 22, 29, 37, 46, 56, . . .}, (12) 

fitted by 

1 - s 

7r~s) 

The growth of the degree is quadratic 

d±± = 1 + 71 (n + 1) /2 (14) 

and the entropy vanishes. 

The few subsidiary entropies we have calculated when a ~ b = c = d also vanish, 
and the degree growth is quadratic. 



9±±{^) = -7—7r (13) 
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8. Example 2: Q4 

We have analysed the so-caUed Qa: lattice equation J2l ^1 ■ The defining relation 
is given by: 

^ ((y[ni,n2] - b) (y[„i,„,+i] -b) -d) ((?/[„i+i,„2] - b) (?/[„i+i,„2+i] - ^) - c^) (15) 
+ B ((y[„i,„2] - a) (y[„i+i,„2] - a) - e) ((?/[„i,„2+i] - a) {y[ni+i.n2+i] - a) - e) = f 
with 

d = (a - &) (c - 6) 

e=(fe-a)(c-a) , , 

f ^ ABC{a-b) 

A{c-b) + B {c-a)^C (a-b) 

We have used our approach for generic values of the parameters, that is to say 
without ifTfill . We find for the fundamental evolutions: 

= {1,3,7,13,21,31,43,57,73,91,111,...} (17) 
fitted with the generating function 

5±±(s) = (18) 
(1 - s) 

The growth of the degree is quadratic, 

4"± = 1 + n(n - 1) (19) 

This indicates integrability of the form lfT5|l ! 

It is interesting to calculate more of the entropies, related to initial conditions 
with a different slope, still for unconstrained parameters. For example 

{d["^2][l]} ^ {1,5,13,25,41,61,85,113,...} (20) 

{d["^2] [2]} =={1,3, 5, 9, 13, 19, 25, 33, 41, 51, 61, 73, 85, 99, 113, . . .} (21) 

Both give zero entropy and quadratic growth, as for the fundamental values. 

9. Example 3: Discrete Sine-Gordon 

A multilinear defining relation for the discrete Sine-Gordon equation can be found 
in [mini: 

?/[ni,n2] 2/[iii + l,n2] y[ni,n2 + l] ^[ni+l^na + l] 

-a (y[«i,n2] y[„i+i,„2+i] - 2/[ni+i,«2] y[ni,«2+i]) -1 = (22) 
For this lattice equation, we find 

{d^"i} = {1,3, 7, 13, 21, 31, 43, 57, 73, 91, 111, . . .} (23) 

as in the previous case (the same quadratic growth, vanishing e±±). 
The subsequent calculation of e[i 2][l] and e[i 2][2] yields 

{^S] [1]} = {!' 4, 11, 21, 34, 51, 71, 94, 121, 151 . . .} (24) 
[2]} ^ {!' 3, 4, 8, 11, 16, 21, 28, 34, 43, 51, 61, 71 . . .} (25) 
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fitted respectively by 

^ ^ (s^ + ,s + 1) (1 — s) 

Which mean vanishing entropies and quadratic growth. 
10. Example 4: Non-isotropic model 

There are cases where the various directions of evolution are not equivalent. The 
entropies e±± are not all equal. Take the simple defining relation (see also |l(i[ll7j 'l: 

y[ni,n2 + l]2/[ni,ri2]y[«i + l.n2] + 2/[rii ."2 + 1] + I.n2 + 1] + Vlrn + Lrii] = (28) 

The sequences of degrees for the fundamental evolutions differ: 

= {1,3, 7, 17, 41, 99, 239, . . .} (29) 
{d[;)^]} = {1,2,4,7,14,28,56,...} (30) 
{d["^_j} ^ {1, 2, 5, 10, 20, 40, 80, . . .} (31) 
{dpi]} = {1,2,4,8,16,32,64,...} (32) 
They fit with the generating functions 

1 + ^ (l-^s)(l + ,s + .^) 

^[^+1 l-2s-s^' ■'^[++1 = T^s ' (^^^ 

1 + 1 



so that 



but 



log(2) (35) 



e_+ = log(2.414...) (36) 

It does not seem excluded a priori to have vanishing entropy in some direction and 
non- vanishing entropy in some other direction, but we have not exhibited any explicit 
example of that yet. 

11. Conclusion and perspectives 

The definitions presented here extend to all dimensions {D > 2), and apply to non- 
autonomous equations, as well as multicomponent systems, provided the evolutions 
are rational. 

We may use a defining relation which is not multilinear if we do not insist on 
having all of the 2^ evolutions described in section (5) . For 13 = 2 we may for example 
accept a defining relation which is of higher degree in j/[„i+i.„2] and y[ni,n2+i]- The 
price to pay is to consider only initial conditions with A1A2 < 0. 

We may also consider defining relations extending over more than one elementary 
cell: this means considering equations of higher order |18j . 
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We have done a number of explicit calculations of the fundamental (and 
subsidiary) entropies, beyond the ones presented here. All the lattice equations which 
are known to be integrable have vanishing entropies. These results, and what is already 
known for maps, suggest that the vanishing of entropy, that is to say the drastic drop 
of the complexity of the evolution, is indeed always a sign of integrability. 

Of course the drop of the degree is related to the singularity content of the 
evolution, as it was for maps. A systematic singularity /factorization analysis is one 
of the tracks to follow. In this spirit, one can undertake a description of all lattice 
equations of a given degree, and for a given dimension. This is the purpose of [T6l[T7| . 



One should also determine which of the properties are canonical, that is to say 
independent of any change of coordinates one may perform on the variables. For maps 
this was just invariance by birational changes of coordinates on a finite dimensional 
projective space. Here the situation is made much more intricate by the infinite 
number of dimensions of the space of initial conditions. 

All the entropies we have calculated explicitly are the logarithm of an algebraic 
integer. There is a conjecture that this is always the case for maps [3 EH) we are 
lead to the same conjecture here. 

Another line of research touches upon arithmetic: a link has been established 
between the algebraic entropy and the growth of the height of iterates, when the 

parameters and the values of y[„j are rational numbers ^niEOlIS]- This applies 

here as well, and will be the subject of further analysis. 
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